DESY 03-111 
FERMILAB-Pub-03/314-T 



Slepton Production 
at e^e~ and e~e~ Linear Colliders 

A. Freitas\ a. von Manteuffel^ and P. M. Zerwas^ 

^ Fermi National Accelerator Laboratory, Batavia, IL 60510-500, USA 
2 DESY Theorie, Notkestr. 85, D-22603 Hamburg, Germany 



Abstract 

High-precision analyses are presented for the production of scalar sleptons, selec- 
trons and smuons in supersymmetric theories, at future e~^e~ and e~e~ linear colliders. 
Threshold production can be exploited for measurements of the selectron and smuon 
masses, an essential ingredient for the reconstruction of the fundamental supersymmet- 
ric theory at high scales. The production of selectrons in the continuum will allow us 
to determine the Yukawa couplings in the selectron sector, scrutinizing the identity of 
the Yukawa and gauge couplings, which is a basic consequence of supersymmetry. The 
theoretical predictions are elaborated at the one-loop level in the continuum, while at 
threshold non-zero width effects and Sommerfeld rescattering corrections are included. 
The phenomenological analyses are performed for e^e~ and e~e~ linear colliders with 
energy up to about 1 TeV and with high integrated luminosity up to 1 ab^^ to cover 
the individual slepton channels separately with high precision. 



1 Introduction 



Supersymmetry [1,2] provides us with a stable bridge [3] between the electroweak scale of 
~ 10^ GeV where laboratory experiments in particle physics are performed, and the Grand 
Unification / Planck scale of ~ 10^^ / 10^^ GcV where all phenomena observed at low en- 
ergies are expected to be rooted in a fundamental theory including gravity. Bridging more 
than fourteen orders of magnitude requires a base of high precision experiments from which 
the extrapolation to the Planck scale can be carried out in a solid way. Such a program 
has already been pursued very successfully for the three gauge couplings which appear to 
unify at the high scale [4]. A parallel program should be carried out in supersymmetric 
theories for the other fundamental parameters [5], including the parameters of soft super- 
symmetry breaking, which may be transferred from a hidden sector near the Planck scale 
by gravitational interactions to our visible world. 

A solid base for these extrapolations can be built by experiments at high-energy e'^e" 
and e~e~ colliders [6-9] which, if operated with high luminosity, will enable us to map 
out a comprehensive and precise picture of the supersymmetric sector at the electroweak 
scale. After the chargino and neutralino sectors [10, 11] have been explored earlier, we 
will concentrate in this analysis on the charged scalar lepton sector of the first and second 
generation, in which mixing phenomena are expected to be strongly suppressed^. [The third 
generation and the neutral sector will be summarized in two later addenda while the colored 
sector will be analyzed in a separate report.] We have elaborated the processes 

e+e-^fitK [^ = L,R] (1) 

and 

e^e~ — > ete~ 

[^,J=L,R] (2) 

e e ^ 

at the level of one-loop accuracy. At the thresholds we have calculated the production 
cross-sections for off-shell particles including the non-zero width effects and the Goulombic 
Sommerfeld rescattering corrections, while in the continuum the supersymmetric one-loop 
corrections have been calculated for on-shell slepton production. 

The threshold production of smuons, mediated by s-channel photon and Z-boson ex- 
changes, proceeds through P-waves, giving rise to the moderately steep [3'^ behavior of the 
cross-sections in the velocity (3 — {1 — 4m? /s)^/^ of the smuons. The accuracy that can 
be reached in measurements of the masses m^^ ^ through threshold scans, is nevertheless 
competitive with the accuracy achieved in the continuum by reconstructing the particles 
through decay products in the final states. Non-diagonal and diagonal pairs of selectrons 
however can be excited in S-waves in e+e" and e~e~ collisions, mediated by t-channel neu- 
tralino exchange, and they give rise to the hnear /3 dependence of the cross-sections near the 
thresholds [13]. This steep onset of the excitation curves allows us to measure the selectron 
masses with unrivaled precision. 

^For a summary of earlier work on this subject see Refs. [7, 12]. 
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Selectron production in the continuum is strongly affected by the electron-selectron- 
gaugino Yukawa coupUngs and, as resuh, they can be determined very precisely by measuring 
the cross-sections for the production processes. In this way the identity of the Yukawa 
couplings {g) with the gauge couplings {g), g — g, a basic consequence of supersymmetry, 
can be thoroughly investigated with high precision, as explored first in Ref . [14] . 

The threshold analyses in this report adopt techniques outlined earlier in Ref. [15]. The 
one-loop calculations in the continuum are performed in the dimensional reduction scheme 
(DRED) for regularization and with on-shell renormalization of masses and couplings. This 
program must be carried out consistently for the slepton sector and the neutralino sector. In 
the loop corrections all sectors of the electroweak supersymmetric model contribute, which 
in general do not decouple for large supersymmetry breaking masses [16,17]. A remarkable 
feature is the appearance of anomalous threshold singularities [18], which show up as discon- 
tinuities in the cross-section as a function of the center-of-mass energy. They are induced by 
specific mass patterns of the particles in the loops [19], which are generally expected to be 
realized in supersymmetric models but which are atypical for Standard Model calculations. 

The phenomenological analyses are performed in the Minimal Supersymmetric Standard 
Model (MSSM), based on the parameters of the mSUGRA Snowmass Point SPSla [20]. 
They include effects of initial-state beamstrahlung radiation as well as the decays of the 
sleptons. All contributions are taken into account that lead to the same final state. They 
have been elaborated at the level typical for phenomenological simulations of processes at 
an e~^e~ linear collider in the TeV range. Besides Standard Model backgrounds the most 
important background channels inside SUSY are taken into account explicitly. The dominant 
standard backgrounds, in particular from W^W~ , ZZ and Z7 production, are eliminated 
a priori by proper cuts adopted from previous experimental studies [21]. At the level of 
precision required here, it is also necessary to include sub-dominant contributions from off- 
shell production of gauge bosons and SUSY particles. 

The final picture is quite exciting: Selectron masses can be determined at an accuracy of 
50 MeV, i.e. in the per-mille range, while the masses of the less frequently produced smuons 
are still accessible at the per-cent level. The same level of accuracy can also be reahzed 
in measurements of the Yukawa couplings of the selectron sector, thus allowing for a high- 
precision comparison with the corresponding gauge couplings. In summa: A high-resolution 
picture of the charged slepton sector in the first and second generation can be drawn by 
experiments at prospective e+e~ and e~e~ linear colliders. 

The report is organized as follows. In Section 2 we summarize the main features of 
slepton production and decay in e^e~ and e~e~ collisions at the Born level. Section 3 
presents the predictions for smuon and selectron production at threshold in detail, leading 
us to the aforementioned accuracies expected from selectron and smuon mass measurements 
in threshold scans. In Section 4 slepton pair production in the continuum is described, and 
exploited finally for measurements of the Yukawa couplings in the selectron sector. Partial 
results had been presented earlier in Ref. [22], while additional technical details can be 
found in Ref. [23]. Spectrum and properties of supersymmetric particles in the reference 
point SPSla, relevant for the present study, are summarized in the Appendix for the sake of 
completeness and the reader's convenience. 
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2 Basics of Smuon and Selectron Production and Decay 



2.1 Notation and Conventions 

In this report we restrict ourselves to the Minimal Supersymmetric Standard Model (MSSM) 
as a well-defined framework. Since the muon and electron masses are very small, the mixing 
among L- and R-smuon and -selectron states, partners of the left- and right-chiral leptons, 
can be neglected and the mass eigen-states correspond to the L,R eigen-states. 

In the other sectors of the MSSM, mixing needs to be taken into account. The MSSM 
requires two Higgs doublets and iJ^i which both acquire non-zero vacuum expectation 
values and v^. The fields mix to form the Goldstone and the physical degrees of freedom 
with the mixing angle 

tan/3 = tiu/tid, (3) 

given by the ratio of the vacuum expectation values. 

The charged higgsinos if^^ and the winos mix to form two charginos xf (^ = 1)2), 

while the neutral higgsinos and the gauginos B,W^ form four neutralino mass eigen- 
states x° {i = 1,2,3,4). 

Apart from the electroweak parameters, the spectrum of the charginos and neutralinos is 
described by three mass parameters, the Higgs/higgsino parameter in the superpotential 
and the soft SU(2) and U(l) gaugino parameters, M2 and Mi, respectively. For the charginos 
the mass term reads ^ 

^m,± = -{W-, H-) X j + h.c. (4) 
where W^, H^^ are the Weyl spinors of the charged winos and higgsinos. The mass matrix 

VV2Mwcos/3 /X J 

can be diagonalized by two unitary matrices U and V according to 

generating the mass eigen-states xf- In the chiral representation, the Dirac spinors xf of 
the charginos are constructed from the Weyl spinors as follows, 

x:=(^^] and xt=(^]- m 

The neutralino mass term in the current eigen-basis is given by 

C^^o = -Ir^Yr + i^-C, r = {B,W\HlHl)\ (8) 
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(a) (b) (c) 

Figure 1: Generic leading-order diagrams for the pair production of smuons and selectrons 
in e'^e~ or e~e~ scattering. 



with the symmetric mass matrix 



/Ml -Mz sw C0 Mz sw \ 

M2 Mic^NCp -Mzc^NSp 

-Mzs^cp Mzc^c/s -/I 

\ Mz Sw sp -Mz cw s/3 -n / 



(9) 



in which the abbreviations = sin /? and cp = cos P have been introduced; sw and cw are 
the sine and cosine of the electroweak mixing angle. The transition to the mass eigen-basis 
is performed by the unitary mixing matrix N, 



N*YN-' = diag(m|o, m|o, m|o, m|o) with Xi = N^^i^l (10) 
The Majorana spinors x° of the physical neutralinos are composed of the Weyl spinors as 

(11) 




Explicit analytical solutions for the mixing matrices can be found in Ref. [10]^ 



2.2 Production Mechanisms 

In supersymmetric theories with R-parity conservation scalar leptons are produced in pairs. 
Since mixing can be neglected, the pairs are built of the current eigen-states with chiral 
index L or R. 

Scalar smuons are produced in diagonal pairs via s-channel photon and Z-boson ex- 
changes in e+e~ collisions, see Tab. 1 and Fig. 1 (a). 

Since the intermediate state is a vector, the hehcities of electron and positron must be 
opposite to each other. By angular momentum conservation the scalar smuons are therefore 
produced in P-wave states. This gives rise to the characteristic (3^ behavior of the excitation 
curves close to threshold, with (3 denoting the velocity of the smuons in the final state. 

^Note that a convention for the chargino mass matrix X different from eq. (5) is used in Ref. [10] . 
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Process 


Exchange i^articles 


Orl")ital \A'a.\T> Tliresliold excitation 


% / 4^1: ^ iJ'RiJ'R / iJ'ti^l 




P-wave oc 


4% ~^ ^R^R 


7,^,X° 


P-wave oc 


^L^R ~^ ^L^L 




P-wave oc 
P-wave oc 


^R^R ~^ ^R 


X 

x° 


S-wave oc /3 
S-wave oc /3 




x° 
x° 


S-wave oc /3 
S-wave oc P 


^R ^R 




P-wave oc P^ 



Table 1 : Classification of smuon and sclcctron production modes in terms of tlie exclianged 
particles, tie orbital angular momentum of tie final state wave function and tie rise of 
the excitation curve near threshold. Specific beam polarization states are required for the 
individual channels. 



The cross-sections for the production of RR and LL smuon pairs by polarized elec- 
tron/positron beams may be written as 



a[ele^ ] = -^p 



1 + 5'L 



s - Ml 



<^[^L Sr ^ i^tfii 



~37 



i + 9i9R 



s-Ml 



(12) 
(13) 



with i = L,R and the left and right-chiral Z couplings 



9l = 



-1 + 24 



w 



2swCw 



9r = — ■ 

Cw 



(14) 



As mentioned before, the polarization combinations with equal helicity of electron and 
positron vanish. The electromagnetic coupling a may conveniently be defined at the en- 
ergy scale y^, incorporating properly the running of the gauge coupling. 

The angular distribution of the smuons follows the familiar sin^ 9 rule so that the new 
particles are produced preferentially perpendicular to the e+e~ beam axis. 

The size of the cross- sect ions for smuons in the mSUGRA Snowmass point SPSla runs 
up to 35 fb and 60 fb for left and right-chiral pairs, cf. Fig. 2. The cross-sections generally 
reach a maximum at s ~ 10 m?, while for asymptotically large energies they scale as 1/s. 

Scalar electrons can be produced, besides the standard photon and Z-boson s-channel 
exchanges, via neutralino x° [j = 1 . . . 4] exchanges in the t-channel, cf. Tab. 1 and Fig. 1 (b,c). 
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thereby generating, in addition to the diagonal, also non-diagonal L/R pairs. Moreover, di- 
agonal and non-diagonal selectron pairs can be generated by t-channel neutralino exchanges 
in e~e~ colhsions, see Fig. 1 (c). 

In contrast to the vectorial s-channel amplitudes, the t-channel neutralino exchanges 
allow for S-wave production at the thresholds. To reduce the total angular momentum to 
zero, electron and positron beams are required with equal helicities. If the helicities are 
opposite, standard vectorial exchanges give rise to the familiar P-wave states, cf. Tab. 1. 

The S-wave production processes are particularly appealing for selectron mass measure- 
ments in threshold scans due to the steep onset of the cross-sections oc (3. In e^e~ collisions 
only mixed selectron pairs, clCr, can be produced in an S-wave, while S-wave production 
of diagonal selectron pairs, crCr and clCl, is possible in the e~e~ mode. Moreover, e~e~ 
collisions provide a nearly background-free environment for selectron studies. 

The Born formulae for selectron production by polarized beams read 



4 4 



s-Ml 



j=i k=i 



(15) 



s-Ml 



f [^ = L/R,-^ = R/L], 



6 



(7[eL e. 



R 



with 



^R 



~37 



1 + 9-i 



s - Ml 



4 4 



-'^Lj 



[i = L/R,-i = R/L], (16) 

(17) 



(T[e+e 



IGttq; 



i=i fc=i 

4 4 

-^Lj ^Rj ^Lfe -'^Rfe /i^*", 



= L/R], 



i=i fe=i 



In 



-2(3 + Aj- In 



A,.-/3 



2 ^x°^xg 
s Aj- ± Afe 



In 



4/3 



Ak-P Aj - /3, 
j 

3^k 



Aj — -{ml. — m|q) — 1 and (3 = ■\J~i--Arr^Js, 



(18) 
(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



where for the case of diagonal selectron pairs, brBr and eLeL, 

2 
-I 
s 

while for mixed pairs 

A.— ^«+<-2m|o)-l and /3 = -J {s - m\ - mfj^ - 4mlmir. (25) 
The matrix 

Xij = [(cw + giSw)Nji + (sw - 5'iCw)-^j2] /V2 (26) 

accounts for the neutralino mixing with N being the neutralino mixing matrix, see (10). 

Since the higgsino components of the neutralino states couple with the small electron 
mass to the electron-selectron system, the exchange mechanism automatically projects on 
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the gaugino components of the neutrahno wave functions. The exchange of relatively light 
neutralinos with dominant gaugino components in the t-channel leads in general to large 
production cross-sections. Only if the neutralinos mass is nearly zero, the contributions 
associated with a Majorana mass insertion in the t-channel of amplitudes with zero total 
spin is suppressed. In the reference point SPSla, however, the masses of the gaugino-like 
neutralinos are sufficiently large, being of the order of the slepton masses, to generate large 
cross-sections in all selectron channels. 

A typical set of selectron production cross-sections is shown in Fig. 3 (a/b) for e+e^ and 
e^e" collisions, respectively. As expected, the t-channel neutralino exchange enhances the 
cross-sections considerably. The cross-sections exceed smuon production by about an order 
of magnitude. 

The angular sparticle distributions for the various production processes and polarized 
beams read 



da ^ - 



a 



4s 



s-Ml 



j=l k=l 

i=i L 



sin^^ 



s 
I0? 



[Aj -/?cos^] [Afc -/?cos^] 
sin^^ 



(27) 



Ml 



Aj - Pcose 



[i = L/R,-i = R/L], 



da 
df)' 

da 
dQ 

da 
dQ 

da 

dn 

da 



\e7 eT 



^R 



^R 



^R 



a 



4s 
16a2 



[i = L/R,-i = R/L], (28) 



myomc,o / s 



s 

da 
dQ 

16a2 



S S ^" - /J COS »] [ A. - ,J cos e] ' 



4 4 



^R ] ' 



4a2 
s 



j=l k=l 

A A 



4A,Afem^om^c/s [i = L/R], 

[A2 _ ^2 eos2 e] [A2 - /32 cos2 e] (30) 



sin^^ 



E E ^R, ^L. Xk. [A^._^,„3^] [a,-/? cose] 



(31) 



with 9 being the angle between the incoming e~ and the outgoing e~ particles and Aj defined 
in eqs. (24), (25), respectively. Near the thresholds the angular sparticle distributions are 
oc sm^ 9 for P-waves while S-wave distributions are isotropic. With rising center-of-mass 
energy, the t-channel neutralino exchange however accumulates the selectrons in the forward 
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Figure 3: Born cross-sections for selectron pair production in unpolarized e^e (a) and 
e~e^ (b) collisions. 
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and backward directions as the exchange amphtudes peak near cos^ ±1: 



dcos^^ ' ^ [Aj- -/3cos^] [Afc -/3cose] 1-cos^' 

da _,_ _ __|_ ^ ■^—v 1 «>"^ij 1 



dcos^' ^ ^ [A,- -/3cos^l [Afc -/?cos^l (l-cos^)^ 

2.3 Decay Mechanisms 

The R-sleptons /xr and cr are expected to decay predominantly into the hghtest neutrahno 
if the latter has a dominant bino component: — > l^Xi- 

The dominant decay modes of the L-sleptons /jl and bl in the SPSla scenario are also 
expected to be decays to the lightest neutralino. However, additional heavy neutralino 
cascade decays and decays to charginos generate more complicated final states [24]. The 
tree-level decay widths for these two-particle decays are given by 

r[C ^ X?] = « m,-^ ( 1 - ^ j - L/R, J = 1 ■ • ■ 4], (34) 



where X denotes the matrix defined in eq. (26), while U is the chargino mixing matrix 
defined in eq. (6). 

Masses, widths and branching ratios for the reference point SPSla [20,24] are collected 
in Tab. 2. While R-sleptons decay almost exclusively into light neutralinos plus leptons, 
the same decay modes are also dominant for L-sleptons. Due to the fairly large value of 
tan/9 and, as a result, the significant stau mixing, charginos and the x^ neutralinos 
decay primarily to r final states, so that their experimental analysis is more demanding. As 
significant rates are predicted for the decay modes of the sleptons into Xii we will focus the 
subsequent phenomenological analyses to these exceedingly clear channels: the final states 
are oppositely charged leptons plus missing energy, e+e~ — > -\- ^ . Decays to x^ with 
subsequent r decays can nevertheless be exploited to discriminate between L- and R-sleptons. 



3 Slepton Production at Threshold and 

Mass Measurements 

Smuon pairs are produced in e+e" annihilation near threshold in P- waves as a result of 
angular momentum conservation for spin-1 photon and Z-boson s-channel exchanges. This 
leads to the (3'^ behavior of the excitation curve in the velocity of the produced particles. In 
contrast, t-channel neutralino exchanges can give rise to a steep linear beta dependence of 
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Sparticle 


Mass m [GeV] 
Width r [GeV] 


Decay modes 




m = 142.72 

r = 0.21 


i^^i-xl 100% 




m = 202.32 
r = 0.25 


!Z ^ I- X? 48% 

^ I- Xl 19% 
^uix^ 33% 




m = 96.18 






m = 176.62 
r = 0.020 


X^^efeT 6% 
^/x^/x^ 6% 
^ 88% 
^ggx? 0.1% 


xt 


m = 176.06 

r = 0.014 


X^^f+i/, 100% 



Table 2: Masses, widths and branching ratios of smuons, selectrons and the hght neutrahno 
and chargino states for the reference points SPSla [20,24]. 



the excitation curves for selectrons in e'^e" and e~e~ coUisions, characteristic for states with 
zero total angular momentum. 

These rules are valid at the Born level but they are modified by the non-zero widths of the 
produced resonances and by Sommerfeld rescattering effects generated by Coulombic photon 
exchange between the slowly moving final-state particles [15]. While the non-zero widths 
smear out the onset of the threshold excitation curves, Coulombic photon exchange enhances 
the cross-section near threshold. For on-shell particle production, the Coulomb correction 
factor is singular, oc (3~^, so that the excitation curves are enhanced to for P-waves and 
they jump to non-zero values for S-waves. For the production of unstable particles, this 
singular behavior is alleviated by the off-shellness and finite width effects. 

Moreover, studying off-shell production of sleptons, the calculation has to be performed 
for the final states after the decays of the resonances. Restricting ourselves to the simplest 
neutralino Xi and X2 decay modes, the processes 

e+e-^//Vx?X? (36) 



and 

e+e-^e+e-x?,2X?,2 (37) 
e~e~ e~e~ x?,2 X?,2 (38) 

must be analyzed including all channels which give rise to these final states. 
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(a) Double resonance diagram (b) Single resonance diagrams 




Figure 4: The doubly and singly resonant contributions to the process e^e — >■ /x"'"//. XiXi- 

The decay of L-sleptons into the next-to-hghtest neutrahno X2 ^'ith the subsequent decay 
X2 'T'^'T^Xi (cf- Tab. 2) can be used to distinguish them from R-sleptons, which predomi- 
nantly decay into the lightest neutralino Xi: i-^- ^ Xi- This is of particular importance 
since in most scenarios the R-slcptons are expected to be lighter than the corresponding L- 
sleptons, so that the L-sleptons are produced on top of a huge background of R-sleptons. 
For scenarios with tan/5 > 10, the X2 mainly decays into a r pair and the lightest neutralino, 
so that the production of an L-slepton is signaled by the appearance of additional r jets. 

The small cross-section for L-smuon production together with the branching ratio for 
fl^ jji^ X2 results in expected event rates that are too low to perform a measurement of 
the threshold excitation curve. Therefore, only R-smuons, but selectrons of both L and R 
type will be analyzed in detail. 

3.1 OfF-shell Slepton Production 

The leading contribution to the ji^ j.!^ x\ Xi final state, and to electron final states corre- 
spondingly, is generated by the double resonance diagram shown in Fig. 4 (a). For invariant 
masses near the smuon mass, the smuon propagators must be replaced by the Breit- 
Wigner form, which explicitly includes the non-zero width F^ of the resonance state. This 
is achieved by substituting the complex parameter 

Ml^ml- im^r^ (39) 

for the smuon mass. To keep the amplitude gauge invariant, the double-resonance diagram 
of Fig. 4 (a) must be supplemented by the single-resonance diagrams of Fig. 4 (b). 

3.2 Coulombic Sommerfeld Correction 

The Coulomb interaction due to photon exchange between slowly moving charged particles in 
Fig. 5 gives rise to large corrections to the excitation curve near threshold. For stable particles 
the cross-section is modified universally by the singular coefficient cTBom («7''/2/3)crBorn at 
leading order. This Sommerfeld correction [25] removes one power of the velocity /3 off the 
threshold suppression. For the production of off-shell particles the singularity is screened [26] 
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Figure 5: Coulomb correction to smuon production, 



X° e+e n+n XiXi- 



and the remaining enhancement depends on the orbital angular momentum /. The Coulomb 
correction is associated with the leading term in the f3 expansion of the photon exchange 
diagram Fig. 5. Exploiting the fact that terms with the loop momentum in the numerator 
generate terms proportional to /3 in the diagrammatic analysis, it follows that 



^oflf-shell 
'^Coul 



— CBorn 7^ Co 3^6 
ZTT 



( 



-2Ml 



\ 2p^p- — p^ — pz 



(40) 



for the complex pole masses M]_ — m\ — im±V± and the momenta p± of the slepton P, 
I = e^jl. The leading part in !3 of the scalar triangle function Co can be evaluated according 
to Ref. [27]. The last factor, taken to the Z-th power, is of kinematical origin, incorporating 
the effect of the angular momentum / of the wave function. 

After carrying out the expansion for smuon and selectron P-wave production and for 
selectron S-wave production, the leading order can be written in the form 



CCoul — CBorn 



1 arctan 



[Z = 0,1], 



with the generalized velocities 

- 



--^(s -Ml- Miy - AMlMl, 
' s —p\— pVf- — 4p^p?. 



(41) 

(42) 
(43) 



The off-shellness damps the singularities as illustrated in Fig. 6 for the S- and P-waves. 



3.3 Final-state Analysis 

The final states in the general processes e+e" — > + If and e"'"e~ /+/^ r+r^ + Ifi can 
be generated by a large variety of background processes in addition to the signal slepton 
channels. Within the SUSY sector itself, pair production of charginos and neutralinos with 
subsequent (cascade) decays feed the final state V^l" + ^. ZZ and Zho/ZHo intermediate 
states with one particle decaying to lepton pairs, the other to a pair of x? will also contribute 
to this class of final states. The final state with an additional tau pair is characteristic 
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270 280 290 300 

^/i [GeV] 

Figure 6: Correction factor Acorr due to Coulomb rescattering relative to the Born cross- 
section for on-shell and off-shell S- and P-wave production. 



for the production of a L-slepton together with a R-slepton. SUSY backgrounds to this 
signature arise from neutrahno production or from stau f production with the decay chain 
e+e" f+f- T+T- X2 Xi T+T- e+e" Xi Xi- 

Moreover, pure Standard Model processes, like the production of gauge boson pairs, 
W^W~ , ZZ and Z'-/*, leading to the final state also have to be taken into account. 

They are generically large and need to be reduced by appropriate cuts [28]. The background 
from resonant Z production can easily be reduced by cutting on the invariant mass of the 
lepton pair or the invisible recoil momentum around the Z-pole. Contributions from WW 
pair production have a characteristic angular distribution of the final state leptons. Because 
of the spin correlations and the boost factor, the leptons tend to be aligned back to back and 
along the beam direction. Therefore this background can be reduced effectively by rejecting 
signatures with back-to-back leptons. The explicit values for the cuts are summarized in 
Tab. 3. 

Triple gauge boson production, W^W~Z and W^W'"-^* , contributes to the final state 
e"^e~ l^l^ T^T^ + Ip . The total cross-section for these processes is well below 1 fb [29] 
and can be reduced further by applying cuts on the invariant di-lepton mass. 

The dominant supersymmetric backgrounds involve decay cascades of neutralinos and 
charginos that, for example following the decay chain 



e^e-^xlx] 



(44) 

i^i-xl 
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Condition 


Vaiial)le 


Aeee])ted range 


ntycCL icpionb ui loiwcuu./ 
backward region from Bhabha/ 
M0ller scattering 


1 OT^'J" j^Ti "i~vi^ 1 T* n rr 1 £i 


1 cob t/1 1 <^ u.yo 


rvejcCL bOiL ibp Lulib / je Lb lioiii 
radiative photon splitting and 
7-7 background 


lcn~\T/~\n / TOT ctt r ni 1 
IcpiUll/Jcl Cllci^y J-Jl 


-C/1 ^ o vje V 


Reject missing momentum in 
forward/backward region from 
particles lost in the beam pipe 


missing momentum polar 
angle 9^. 

*— ' i^miss 


Icos^p- . 1 < 0.90 

1 /'miss 1 


Angular separation of two 
leptons 


angle 0i+i- between leptons 


|1 - cos 01+1-1 > 0.002 


Angular separation of two quark 

or tau jets 


angle 0jj between jets 


|1 - cos0jj| > 0.015 


Cut on Z decaying into 


di-lepton invariant mass 


|mi+i- - Mzl > 10 GeV 


lepton pair 


mi+i- 




Cut on invisibly decaying Z 


invariant recoil mass mrecoii 


Imrecoii - Mzl > 15 GeV 


Reject back-to-back 
leptons from W pairs 


angle 0i+i- between leptons 


1 cos 01+1- 1 < 0.7 



Table 3: General cuts to reduce large Standard Model backgrounds and to account for the 
detector geometry and resolution. 



with j > 1, generate -|- final states. Since the lepton pair originates only from a 
single neutralino decay, these backgrounds give rise to increased missing energy and lower 
lepton-pair invariant mass compared to the signal, and they can effectively be reduced by 
cuts on these two variables [15]. Near threshold an alternative method for reducing the 
backgrounds can be applied, based on the fact that the energy of the leptons originating 
from a two-body decay is defined sharply in this kinematical configuration. Thus by selecting 
leptons with energies in a band AE ?a 10 GeV around the nominal threshold energy -E;,thr ~ 
(m~ — m'iq)/{2mj) greatly suppresses both SM and SUSY backgrounds. This second cut 

choice is applied in the following examples. 

The signal-to-background ratio can further be enhanced by using beam polarization. The 
optimal polarization choices for the different production processes are listed in Tab. 1. As 
evident from the table, polarization of both the electron and positron beams can help to 
discriminate between the slepton chiralities (see e.g. [30]). In the following 80% polarization 
for the electrons and 50% polarization for the positrons is assumed. Without positron 
polarization the signal-to-background ratio would in general be reduced by a factor of 1.5. 

For a realisitic description, it is necessary to include initial-state radiation (ISR) and 
beamstrahlung effects. The leading logarithmic ISR contributions are included using the 
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^ [GeV] 



Figure 7: The excitation curves for cr pair production over Standard Model and super- 
symmetric backgrounds for e+e~ anniliilation (left) and e~e~ scattering (right). The signal 
contribution with non-zero widths and Coulomb rescattering is compared with the case of 
zero width and no higher order corrections. The signal is enhanced with beam polarization 
as indicated, where (+) corresponds to right-handed and (— ) to left-handed polarization. 



structure-function method [31], while effective beamstrahlung parametrizations, taken for 
the Tesla design for definiteness, are adopted from the program Circe [32]. 

Including the SUSY and SM backgrounds, the excitation curves, after the beamstrahlung 
and ISR is switched on and the cuts defined before are applied, are displayed in Fig. 7 for 
two characteristic examples, e^e^ pair production in e^e~ collisions as a P-wave process, 
and Cj^Cj^ pair production in e~e~ collisions as a typical S-wave process. Separately shown 
are the zero-width Born prediction, the background contributions and the final prediction 
including non-zero width and rescattering effects, with the backgrounds added on. 

The results expected from these simulations for the mass measurements are presented 
in Tab. 4. They are based on data simulated at five equidistant points in a center-of-mass 
energy range of 10 GeV in the threshold regions for ftn pair production, and diagonal and 
non-diagonal cr and cl production. For the e~^e~ mode a total luminosity of 50 fb~^ for each 
threshold scan is assumed, corresponding to 10 fb~^ per scan point. In the e~e~ mode the 
anti-pinch effect leads to a somewhat reduced machine luminosity. Therefore it is presumed 
that a total of 5 fb~^ is available for each scan measurement, corresponding to 1 fb~^ per 
scan point. For the reconstruction of the mass a binned likelihood method is employed, 
using four free parameters in the fit: the slepton mass and width, a constant scale factor for 
the absolute normalization of the excitation curve and a constant background leveP. The 

■^Since the remaining backgrounds after cuts are flat, they can efiiectively be approximated by a constant. 
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Process 


Fitted values for slepton mass m and width F 


e+e- ^ (e+gj,) ^ e+e" + ^ 


mg^ = 142.8l°:2i GeV T,^ = ISO^^oo MeV 


e-e" ^ (e^eR) ^ e'e" + ^ 


rrie^ = 142.70l|J:[^^^ GeV Fe,, = 200tlo MeV 


e+e" (e^e^) e+e" r+r" + ^ 


rrie^ = 202.21°!^ GeV Fg^ = 240i^g MeV 


e"e- — > (e^e^) — > e"e" tttt + ^ 


= 202.11°:^^ GeV Fg^ = 2401^^^ MeV 


e+e" ^ (/i^/iR) ^ + ^ 


m^, = 142.8l°:i GeV F^, = 35011°° MeV 



Table 4: Expected precision for the determination of slepton masses and widths from 
threshold scans in e+e" and e~e~ scattering. The reconstructed values are obtained from a 
four-parameter ht as outlined in the text. 

last two parameters render the mass fit independent on details of other SUSY sectors, in 
particular the masses and mixings of the heavier neutralinos that are not accessible in the 
slepton decays. 

Evidently, S-wave production in e-e" collisions provides us with mass measurements 
of 50 MeV, i.e. a relative error of less than 1 per-mille. This will presumably be the highest 
accuracy that can ever be reached for sfermion mass measurements in the supersymmetric 
particle sector. 

4 Slepton Production in the Continuum and 
Determination of Yukawa Couplings 

The motivation for high-precision analyses of smuon and selectron production in the contin- 
uum is twofold, different though for the two species. 

1. ) Smuon pair production in the continuum serves as a rich source of particles which 
after subsequent decays to muons and neutralinos allows us to determine the smuon and 
neutralino masses. The smuon mass measurement in the continuum is competitive with the 
accuracy expected from threshold scans. The R-smuon decay to the lightest neutrahno will 
be the gold-plated process (besides the analogous R-selectron decay) for measuring the mass 
of the lightest neutralino, which is a key particle in cosmology. Smuon pair production leads 
to clean yU+yU^ + ^ final states that can easily be identified and controlled experimentally. 
The measurement of the cross-sections therefore provides a valuable instrument for testing 
supersymmetry dynamics at the quantum level. 

2. ) Since the precision of threshold scans in mass measurements of selectrons cannot be 
rivaled, the central target of selectron pair production in the continuum, besides the neu- 
tralino mass measurement [21,33], is the analysis of the sclcctron-clcctron-neutralino Yukawa 
couplings in the SU(2) and U(l) sectors [16]. They are predicted to be equal to the corre- 
sponding gauge couplings in supersymmetric theories, even if the supersymmetry breaking 
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is included by soft terms in the Lagrangian, leaving the system theoretically self-consistent. 
The relevant mechanism involves the t-channel exchange of neutralinos. Knowledge of the 
neutralino masses and mixing parameters is therefore required before high-sensitivity tests 
can be carried out. Thus this method is one of the components in a complex experimental 
program including, in addition to the analysis of selectron pair production, also pair produc- 
tion of charginos and neutralinos that in turn are (partly) mediated by neutral and charged 
slepton t-channel exchanges. The synopsis of all these channels will finally provide us with 
a comprehensive and detailed picture of the entire Yukawa sector in a model-independent 
form. 

We will separate in this report the description of the theoretical techniques necessary for 

controlling the higher-order corrections, from the application to smuon and selectron pair 
production and their phenomenological evaluation, with emphasis on the analysis of Yukawa 
couplings. 

4.1 Renormalization of the MSSM 

At the one-loop level, which we work out in this report for slepton production, dimensional 
reduction (DRED) provides us with a valid regularization scheme including chiral currents^. 
By reducing the kinematics in the propagation of particles to D < A dimensions but leaving 
the number of field components unchanged, supersymmetry is preserved in the higher-order 
amplitudes, and so is gauge invariance [34]. 

Multiplicative renormalization of masses, couplings and fields can therefore be performed 
without introducing additional ad-hoc counter terms to restore the supersymmetry. The 
renormalization factors Z, and equivalently the shifts of variables, that absorb all the ultra- 
violet divergences, will be fixed by on-shell renormalization, i.e. the on-shell definition of the 
physical particle masses, the on-shell definition of the electromagnetic gauge couplings in the 
trilinear lepton-lepton-photon vertex, and normalization of the on-shell renormalized fields 
to unity. As a consequence of supersymmetry and gauge symmetry, the renormalization of 
all other quantities, Yukawa couplings, quartic couplings etc., induces calculable additional 
shifts. This program can be carried out consistently in theories including soft supersymmetry 
breaking terms^. 

Characteristic classes of higher-order diagrams for propagators and vertices are depicted 
in Fig. 8. Additional box diagrams, cf. Fig. 9, finally conclude the set of elements contributing 
to the 2-2 transitions^. 

After carrying out the renormalization program in the ultraviolet sector, infrared and 
coUinear divergences associated with the massless photon and lepton fields can be absorbed 
by adding the real photon emission contributions. Fig. 10. Finite results are automatically 

^For tli(^ production of smuons, involving only gangc couplings at tree level, the calculation has been 
repeated independently in dimensional regularization and perfect agreement has been found [23]. 

^More information about the renormalization of the MSSM and a general overview of different renormal- 
ization techniques can be found in Ref. [35]. 

®Most of the analytical results for self-energy operators etc. are too lengthy to be presented in this report; 
therefore computer codes for the calculation of the loop results are made available on the web, cf. the 
concluding remarks in section 5. 
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Figure 8: Generic sets of Feynman diagrams for the virtual self-energy and vertex cor- 
rections to slepton production. Solid, dashed and wiggly lines indicate fermions, (Higgs) 
scalars and vector bosons, respectively, whereas sfermions and gauginos are denoted by dou- 
ble lines and wiggly/solid lines. The selectron-electron-neutralino vertex in the last line only 
contributes to selectron production. 
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Figure 9: Box-type Feynman diagrams for slcpton pair production. The first row (a) 
applies botli to smuon and selectron production, wliile tlie second row (b) only contributes 
to selectron production. 



(a) (b) 




Figure 10: Feynman diagrams for real plioton emission in smuon production (a) and selec- 
tron production (a)-h(b). 



guaranteed by proceeding to experimentally well defined cross-sections, i.e. the total cross- 
sections in the present analysis. 

Due to large number of diagrams involved, the use of computer algebra tools for the 
computation is necessary. The generation of diagrams and amplitudes is performed with 
the package FeynArts [36] . Throughout the calculation, the CKM matrix is taken diagonal 
and mixing between the sfermions of the first two generations is neglected. For the third 
generation sfermions, the mixing between the L- and R-statcs is consistently taken into 
account. A general covariant gauge is used in order to facilitate an additional check of 
the result. Using the program FeynCalc 2.2 [37], the Lorentz and Dirac algebra is evaluated 
and the loop integrals are reduced to a set of fundamental scalar one-loop functions [38]. 
Since the explicit analytical expressions of the virtual loop contributions are generally very 
lengthy, they have been implemented into a computer code that calculates the one-loop- 
corrcctcd cross-sections, using the package LoopTools [39] for the numerical evaluation of 
the basic scalar one-loop functions. 



4.1.1 Gauge sector 

The extension of the Standard Model to a supersymmetric theory in minimal form (MSSM) 
does not introduce new couplings in the gauge/gaugino sector. The gauge sector of the 
MSSM is therefore renormalized in parallel to the Standard Model. Just the self-energies 
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are expanded by the contributions of the supersymmetric fields in a straightforward way. 
We briefly summarize the results, adopting the standard conventions of Ref. [40]. 
The masses of the W and Z gauge bosons are shifted by 

Mlr + SMlj and , , xa#2 



Ml 



Ml 



,W -r — -^Z - -^Z -r SM'l (45) 

Imposing the on-shell renormalization conditions defined earlier, the mass shifts can be 
expressed in terms of the transverse self-energies St, 

SMl = S^W(M^) and 6 Ml = (M|), (46) 

with the self-energies for the gauge boson propagation in the Standard Model expanded by 
supersymmetric particle contributions: 

E^^^^(^2) = E^^^^(fc2)|3^ + E^^^^(fc2)|3^3^ [V, = 7,Z, W]. (47) 

The first term includes the usual Standard Model loop contributions as in the first line of 
Fig. 8, while the second term accounts for the additional loops involving pairs of supersym- 
metric fields, gaugino and sfermion fields, as given in the second line of Fig. 8. 

The SU(2) and U(l) gauge couphngs g and g' can be traced back to the electromagnetic 
coupling e and the electroweak weak mixing angle sw = sin^w e is renormalized as in 
standard QED apart from the removal of ^y-Z mixing, 

1 dE:p{e) sw ^1^(0) 

2 



e ^ (1 -|- 6Ze) e with 



fc2=0 ^ 



M| 



(48) 



Again the gauge-boson self-energies decompose into Standard Model and specific supersym- 
metric contributions as in (47). 

Introducing the electroweak mixing angle in on-shell definition through the W and Z 
masses as = 1 — M^/M|, the renormahzed value is formally related to the bare value by 



w 



with 



6s 



w 



6Ml 
Ml" 



SM^ 
Ml 



Finally, the renormalized left- and right-handed electron fields, 

CL ^ (1 + I^Z^^) CL, eR ^ (1 + i5Z^^) eR, 

are related to the electron self-energies by 



5Z^^^ 



^^'(^e) + ^ [S^'(/) + S^^(/) + 2/me E^«(/)' 
T,'\ml) + ml ^ [e^^(p') + E^^(p2) + E'^ip^) 



(49) 

(50) 

(51) 
(52) 



(53) 



with the decomposition 

S^(j9) = ^ cuLS^^(p2) + ^ curS^^(/) + S^S(p2) with a;L,R = (1 ± 75)/2. 

Apart from the calculation of the (singular) QED corrections, the chiral limit of vanishing 
electron mass can be safely applied, simplifying eqs. (51), (52) to 

(54) 



-^e E;^^k(0), 



"^weak — 



-^e e;^,,(o). 



21 



4.1.2 Sfermion sector 



In the limit of vanishing lepton masses in the first and second generation, the L- and R- 
selectron and smuon fields do not mix and the mass matrices are approximately diagonal. 
The "chiral" L and R states coincide with the mass eigen-states. This remains true in 
higher orders as the sfermion mixing is proportional to the associated lepton mass. The L 
and R fields may therefore be treated independently so that the renormalization follows the 
standard procedure. With 

rnf — >• mj + Smj and — >• (1 + ^SZ^) li [I = e, iJ,;i = L, R], (55) 
we find for the sfermion mass shift 

5m? = 3f?eE[(m?), (56) 
and for the wave-function renormalization 



(57) 



Here S-(A;^) denotes the self-energy for the slepton l^; I — e, ji; i — Since the external 

fields are super partners, the slepton self-energies cannot be separated into a Standard Model 
and a genuinely supersymmetric part. 

Note that the mass shift and wave-function renormalization for L-sneutrinos coincide 
with those of L-selectrons in the chiral limit we consider in this report. 



4.1.3 Chargino and neutralino sector 

The spectrum of two charginos and four neutralinos in the MSSM is described by the three 
mass parameters M2 and Mi, see section 2.1. Apart from other electroweak parameters, 
the system is also affected by the Higgs mixing tan/?. Three chargino/neutralino masses 
are sufficient to fix the mass parameters /i, M2 and Mi. The renormalization of tan/9 is 
performed outside the chargino and ncutrafino sector, as will be discussed in section 4.1.4 

The other three masses and the mixing parameters are then uniquely determined once 
the parameters in the loop corrections are known [41,42]. Following [41], the renormalization 
is performed in the current eigen-basis. 

1.) Starting from the chargino Lagrangian 

>Cch = ^[V'-V5^^ + ^^a^9^V'+] - [V^-^X^+ + V^^XtV^], (58) 
with the current fields 
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the mass matrix X is renormalized by 
X ^ X + 5X with SX 



SM2 V2 5(Mwsin/3) 

V2 5 (Mw cos (3) 5n 

and the current fields are replaced by the normalized mass eigen-fields x^, 



+ 



X 



+ 



4^- 



X 



(60) 



(61) 



Besides the renormalization of the new parameters /i and M2, SX includes the renormaliza- 
tion of tan /? and the W mass discussed earlier. The [infinite] multiplicative renormalization 
of the wave functions is absorbed in the matrices 5Z^, SZ^, and so is the [finite] renormaliza- 
tion of the matrices V, U rotating the current to the mass fields. The renormalized chargino 
Lagrangian and the associated counter terms may be written in the form 



XI 

xt 



5C 



ch 



xt^ xt 



+ 5Z^) + 1- {SZ^'' + ' ) uJK 



- {\5Z^^U*XV^ + \U*XVHZ^ + WdXV"^) uji^ 

- (l^z^Vxtc/^ + ivx^u^sz^* + VSX^U^) 



(62) 



(63) 



XJ 

xt 



The physical x"^ masses can be introduced in (62) after diagonalizing this part of the La- 
grangian by rotation through U, V. The countcrterms 5fi and 6M2 can thereby be adjusted 
such that the propagator matrix develops poles at the on-shell chargino masses iT^^f^- In 

addition, the Z^'^ factors can be uniquely fixed by requiring that the propagator matrix is 
diagonal and that the pole residues are normalized to unity for on-shell momenta. 

2.) The analogous program can be carried out in the neutralino system, though the 
doubling of degrees of freedom renders the analysis more cumbersome. The bilinear part of 
the neutralino Lagrangian in the current eigen-basis is given by 



1 



(64) 



I 

with and Y given in (8) and (9), respectively. In this representation the renormalization 
of the mass matrix Y is defined as 



Y ^Y + 5Y, 



with 



5Y^ 



( SMi -5(MzSwc^) S{Mzs^^sp)\ 

SM2 S{Mzc^Cf3) -5{Mic^sp) 

-6{Mzs^^cp) S{MzCwC(3) -5/1 

\5(MzSwS/3) -5(MzCwS/3) -Si^ / 



(65) 
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while the renormahzation and the rotation from current fields ijj^ to mass eigen-spinors x° 
can be performed as 

^N^ (l + x°- (66) 

The matrix SZ^ absorbs the multiplicative renormahzation of the current fields as well as 
the renormahzation of the rotation matrix. Thus the renormalized neutralino Lagrangian 
can be cast into the form 



X? 



(67) 



together with the counter terms 

1 
2 



{\5Z^^N*YN'^ + \N*YN^Z^ + N*5YN'^) 
{\8Z^^NY^N'^ + |iVFtiV^5Z°* + N5Y^N^) cjr 



(68) 



X? 
Xl 



where the matrix rotates the neutralino mass matrix into diagonal form according to (10). 
The mass m^o of the lightest neutralino that will be under excellent experimental control, 
may be chosen to define the remaining U(l) gaugino mass parameter Mi. The masses of 
the heavier neutralinos are thereafter fixed uniquely by the Higgs/higgsino and gaugino 
parameters and M2,i. Again, the elements of the Z^ wave-function renormahzation matrix 
can be adjusted such that the elements of the neutralino propagator matrix are diagonal 
with unit residues of the mass poles for on-shell momenta. 

In the case of CP conservation, the renormahzation of the Higgs/higgsino parameter fi 
and the gaugino parameters Mi and M2 may be cast in the following form 



5Mo 



5^ 



SMi 



^±1, - m^±M2) 5Re{m-±S±L(m^±) + m-±E±^(m^±) + 2 ^^^^^(mji)} 
+|(m-±/x - m-±M2) Ke{m-±E±^(m2±) + m-±E±^(m^±) + 2 ^tt{ml^)} 
+M2 5Ml + n 5{MIj sin 2/3)] / (/x' - M|), 

'i(m-±M2 - m-±/.) 3?e{m-±E±L(m^±) + m^^ll^^imi^) + 2 S±SL(m2±)} 
+|(m^±M2 - m^±^x) 5?e{m^±E±^(m^±) + m-±E|^(m^±) + 2 E|SL(m^±)} 
+// 5714, + M2 (5(M^ sin 2/5)] / {Ml - fi^), 

^ - fepim-oEOL(m|o) + E?f (m|o)} - N^^ SM2 + 2NisNu Sfx 
+ 2Nii [Ni^ 5(MzSwCOS/3) - 5(MzSwsin/3)] 



(69) 



(70) 



Nil 



(71) 



+ 2Ni2[Ni^ 5(MzCwcos/?) - Ni^ 5{M^c^sml3)] 
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which is in agreement with [41]. Here the following decomposition of the chargino/neutralino 
self-energies has been used: 

E?,(p) = j^CLEf (/) +^u;n^fip')+u;^J^f^ip')+u;n^''ip'). (72) 

The combinations of self-energies in (69)-(71) are equivalent to the counterterms of the 
on-shell chargino and neutralino masses, 

Sm-^ = ^5Re{m-±E±^(m^.) +m-±E±^(m^,) + 2E±f (m^.)} [k = 1,2], (73) 
6m^o= Ke{m^oE?5^(m|o) + E?f(m|o)}. (74) 

Once the two Xi,2 chargino masses and the Xi mass are fixed, the remaining heavier neutralino 
masses are shifted by finite amounts relative to the Born terms [which, by definition, are the 
eigenvalues of the renormalized mass matrix [35,41]], 

^xE - = -M^xl^'t^Kl) + ^IfKl^} + iN*^yN^)kk [k = 2, 3, 4]. (75) 

The cancellation of the divergences between the neutralino self-energies and the mass matrix 
counterterm SY in this expression is a non-trivial check of the method. 



4.1.4 Higgs mixing tan /3 

At tree level, the Higgs mixing parameter tan /? is determined by three soft SUSY breaking 
parameters, the diagonal mass parameters ml and ml as well as the H^-H^ mixing term m|, 
that is connected with the soft parameter Bfj,. They define the bilinear part of the scalar 
potential of the two Higgs doublets, 

Hilin = ml (10^ + yi + |0d I') + ml (10^ + yi + ^^^^ 
+ml (0u0d + PuPd + 0^0J + 0r<^d *) ■ 

The three soft SUSY breaking parameters can be reexpressed in terms of the vacuum expec- 
tation values and v^, and the pseudoscalar mass M^- The tree-level relation tan /3 = v^/vd 
is modified by higher-order corrections to the Higgs potential and needs to be renormalized. 

Various renormalization prescriptions have been proposed in the literature, which however 
do not lead to satisfactory solutions on all accounts [43] . Renormafization prescriptions that 
are derived from the Higgs potential introduce dangerously large corrections to tan (3, which 
effectively invalidate the convergence of the perturbation series. Other methods impose the 
renormalization condition that the Goldstone bosons to not mix with the physical Higgs 
bosons for on-shell momenta [44], but the value of tan/3 in these schemes depends on the 
gauge choice. 
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Alternatively, one may define tan f3 in higher orders by relating it to a specific physical 
process. However, any definition of tan f3 through a physical process is afflicted with technical 
difficulties that are introduced by the particular process. Moreover, the value of tan (3 may be 
extracted from different observables, for example from Higgs decays [43,45], Higgs production 
processes [46], or from the cross-sections for mixed chargino pair production for moderate 
values of tan/? < 10 [10], while for large tan/? the r polarization in f decays provides an 
attractive opportunity [47]. Any such approach does not lead to a unique and universal 
choice for the renormalization of tan /?. 

In the following analysis we adopt the most convenient solution by just subtracting 
the divergent part in the DR scheme from the unrenormalized parameter to define the 
renormalizcd parameter. Without loss of generality, the universal DR counterterm for tan /3 
may be extracted from the mixing self-energy of the Z-boson and the pseudo-scalar Higgs 
boson A^, with the natural choice for the renormalization scale being the A'^-mass Ma: 



with the subscript "div" indicating that only the divergent part of the self-energy is retained. 

Though being process- independent, this definition is not perfect either as the value de- 
pends on the chosen gauge. [By accident it remains independent on the gauge fixing param- 
eter ^ in the gauge to one-loop order]. Of course, the predictions for physical observables 
remain gauge independent as the gauge dependence of tan/^jj^^; is neutralized by equivalent 
terms in the amplitudes themselves. 

4.2 Effective Yukawa Couplings 

In general, quantum corrections are reduced with increasing mass of the virtual particles 
inside the loops, as generally expected by the uncertainty principle and formalized by the 
decoupling theorem [48] . However, in theories with broken symmetries, these corrections may 
grow to large values if the mass splitting in the particle multiplets is large. High mass scales 
in theories with broken symmetries can thus manifest themselves in the radiative corrections 
to precision observables at much lower energies. A classical example for this phenomenon 
is provided by the mass splitting in the top-bottom iso-doublet of the Standard Model [49] 
which strongly affects the p parameter in the ratio of W- and Z-boson masses. 

In theories with broken supersymmetry such a phenomenon arises if the splitting be- 
tween the masses of SM particles and some of the SUSY partners becomes very large [such 
scenarios are realized, for instance, in focus point theories [50]], leading to superoblique cor- 
rections that grow logarithmically with the mass splitting [16, 17]. The splitting affects in 
particular the relation between the gauge couplings and the associated Yukawa couplings. In 
parallel to the bare couplings, the renormalized Yukawa couplings can naturally be defined 
to be the same as the gauge couplings at the renormalization point if supersymmetry is 
broken softly. Accordingly, in the on-shell renormalization scheme the renormalized Yukawa 
couplings are defined to be equal to the renormalized on-shell gauge couplings, so that the 
renormalized Lagrangian manifestly reflects the supersymmetry [modulo the soft-breaking 



tan f3 — > tan ^ + 5 tan /3 



with 
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Figure 11: Quark-squark loop corrections (a) to the U(l) and (b) to the SU(2) neutral 
gaugino self-energies. 

terms]. Nevertheless, the loops will modify in toto the two types of vertices associated with 
the two couphngs differently in physical amplitudes evaluated near the light SUSY scale, 
or equivalently the electroweak scale. In a more intuitive language, the running of the two 
couplings with energy from the high SUSY mass scale down to the low SUSY mass scale (or 
electroweak scale) is different [16,17]. 

These effects are rooted in the self-energies of the gaugino hnes, shown in Fig. 11 (a) and 
(b) for bino and wino lines, respectively. In parallel to the iso-multiplets of the Standard 
Model, the non- decoupling corrections arise from large mass splittings within the supersym- 
metric particle spectrum, for example for exceedingly high squark masses. The wave-function 
renormalization associated with these diagrams can be projected onto effective U(l) and 
SU(2) Yukawa couphngs, g'^g and ^eflf, defined near the light SUSY/electroweak scale Me^; 
in the same way as the wave-function renormalization Z'^ of the photon defines the effective 
electromagnetic coupling e^s = \f7n e. 

The self-energy correction derived from Fig. 11 (a) for the U(l) bino propagator is given, 
in leading logarithmic order of the ratio between the large SUSY and the electroweak scale, 

by 

4,.o,(P) = ^^ylnM^/Mi. (79) 
This may be reinterpreted as a shift of the effective U(l) Yukawa coupling, 

in leading logarithmic approximation, in relation to the bare Yukawa/gauge couplings g'^ — 
(yfQ. The effective U(l) gauge coupling g%^^a^ can be introduced in parallel. It is identical to 
the renormalized gauge coupling g' in the on- shell scheme. 

The ratio of U(l) Yukawa to gauge coupling is therefore given effectively by 
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: Vertex graph for the process e~^e~ fl'^fl~ 
an anomalous threshold for the mass hierar- 



Thus, even in spite of the fundamental identity of the renormahzed Yukawa and gauge cou- 
pUngs in supersymmctric theories in soft SUSY breaking scenarios, we find nevertheless log- 
arithmically enhanced departures from universality in effective couplings at the electroweak 
scale if the mass splitting in the supersymmetric multiplets is large. 

The SU(2) Yukawa and gauge couplings can be treated analogously. From the wino self- 
energy in Fig. 11 (b) the effective SU(2) Yukawa couphng g'^^ can be defined. Together with 
the on-shell renormahzation of the SU(2) gauge coupling g, this leads to 



:;2 ^2 



g = 1 + 3 In M?/M,^^ (83) 

for the ratio of the two effective Yukawa and gauge couplings in leading logarithmic order. 

The logarithmic growth of the one-loop corrections will later be analyzed numerically 
when the production of selectron pairs, involving the Yukawa couplings in t-channel neu- 
tralino exchange amplitudes, will be compared with the production of smuon pairs in detail. 

4.3 Anomalous Thresholds 

The rich pattern of different masses in supersymmetric theories gives rise to anomalous 
threshold singularities [18, 19] in vertex and box graphs [which play no role in general in 
the Standard Model^]. While the vertex graph for e'^e~ — > in Fig. 12 generates a 

normal threshold singularity when the energy y/s passes the threshold for XiXa production, 
an additional anomalous singularity occurs for a special set of mass values m^o < < m^o 
at the kinematical point 

s^s = ^ f84) 

Here, as before, the muon mass has been neglected. The singularity can be traced back to 
a zero value of the denominator function D in the vertex amplitude /, in a configuration 
where all intermediate particles in the loop become on-shell at the same time, 

I- Jd'ql^^l^^l^, D^q'[iq + p,.r-mlo\[iq + p,.-p,-r-mlo\, (85) 



''As an exceptional case, anomalous thresholds do occur in W^W scattering in the Standard Model 
[comment by T. Hahn]. 
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where / is a polynomial function. Even though the singularity is mild and can be integrated 
over, it leaves its trace in a discontinuity of the cross-section. 

The two different types of threshold singularities for the vertex introduced above are 
exemplified in Fig. 13 and Fig. 14. The normal threshold singularity at ^/s^ in Fig. 13 is 
generated by the non-zero onset of the imaginary part of the 2-point function Bq{s, m|(,, m|o), 
corresponding to a kink in the real part. In contrast, the anomalous threshold at in 
Fig. 14 is generated by the Co(s, , m? , m|„, m|o, 0) function, characterized by a kink in 
the imaginary part which leads, if integrated out in a dispersion relation, to a step in the 
real part [18]. 



4.4 Production Cross-Sections 

Mediated by the pure s-channel 7/Z exchange mechanisms, the production of smuon pairs 
is the most basic process of supersymmetric theories at e"'"e~ colliders. The results for the 
production cross-section will therefore be presented for this process first. Subsequently we 
expand the analysis to selectron pair production in e^e~ and e^e" collisions which involve 
also t- channel neutralino exchange mechanisms. They will serve as an excellent instrument 
to measure the electron-selectron-gaugino SU(2) and U(l) Yukawa couplings, which will be 
analyzed at the end of this section. 



4.4.1 Smuon production 

After the renormalized transition amplitude for the process e^e^ /^^/^ is constructed 
following the way outlined in the last section, the experimental parameters must be defined 
properly. We use the on-shell definition for all masses, while the electromagnetic coupling 
a will be evaluated at the scale of the center-of-mass energy Q = y^, so that the large 
logarithmic corrections oc log s/m^ from light fermion loops in the running of ct{Q'^) are 
absorbed into this definition. 

The resulting amplitude is UV finite but still infrared divergent. This divergence is 
removed by adding the contributions from photon radiation in the initial and final states 
to the cross-section. The virtual and real QED corrections form a gauge invariant subset 
separate from the other virtual loop corrections. 

Initial-state QED corrections: Adding to the loop-corrected cross-section the contribution 
of soft photon radiation from the initial lepton lines, the ensuing cross-section factorizes into 
the Born cross-section and a radiation coefficient that depends on the cut-off AE of the soft 
photon energy (defined in the cms frame). 



J ^virt+soft ^ ^ ^ 



TT 



log^|log^-l +-log^-2 + - 



e 



(86) 



This AE dependence is removed if the radiation of hard photons is added to the cross- 
sections. We are still left, however, with the logarithmic enhancement of the cross-section 
from coUinear radiation oc logs /ml. In leading logarithmic order the photon radiation 
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Figure 13: Effect of tlie normal threshold in the finite part of the 2-point function 
i?o(s, m?o, "n^-o) at the branch cut ^/s^ = m^o + m^o = 474.05 GeV [in the SPSla scenario]. 
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Figure 14: Effect of the anomalous threshold defined in eq. (84) in the 3-point function 



Co(s,m?,m?,m|o,m|o,0) at 



Sa = 516.56 GeV [in the SPSla scenario]. 
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effectively just reduces the cms energy available for the final-state particles [51]. This is 
described by the convolution 



dcrLL(s) = / dz rgg^(2Q;, z, s) CrBorn(^s) 



(87) 



of the Born cross-section with the radiator function [52] 



AVI ^ J_ " 1 1+-^" 



27r 



1-z 



z<l-e 



(that can easily be generahzed to higher orders [53]). The variable z denotes the energy 
fraction left to the electron/positron parton after the radiation of the coUinear photon. The 
additional non-collinear photon radiation is treated numerically by applying Monte Carlo 
integration techniques, with fast convergence after the leading logarithmic order is subtracted 
analytically, as outlined above. 

Final-state QED corrections: After adding up the vertex correction plus the final-state 
photon radiation, the cross-sections for soft photons factorizes again in the Born cross-section 
and a radiation function that depends on the photon cut-off energy A£^, 



J virt+soft 



da 



a 



Born 



TT 



log 



4(AE)2 
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-1+ log 



2/3 
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l-/3V'"2^°'l-/3^ 



(89) 



Lio 



2/3 • 
1 + /3, 



Since the smuon mass is large, the velocity (3 — {1 — 4m? /s)^/^ of the smuons in the final 
state stays sufficiently away from unity not to generate coUinear singularities. The A£' 
dependence is neutralized when the hard photon contributions are added and integrated out 
(numerically) to calculate the total cross-section. 

The initial and final-state QED corrections do not interfere in the total cross-section 
as a consequence of CP invariance. However, the amplitudes do in general interfere in the 
calculation of final-state distributions. 

In general it is not possible to divide the virtual loop corrections for slepton pair pro- 
duction into SM-like corrections and genuine supersymmetric corrections. However, for the 
special case of /ir pair production, a gauge-invariant and UV-finite subset of SM-like loop 
contributions can be defined, including all diagrams where a SM fermion, gauge boson or 
the hghtest Higgs boson is attached to the tree-level graphs, and taking the mixing angle 
a of the CP-even Higgs bosons to be a = (3 — 7r/2. The remaining loop contributions can 
then be interpreted as the genuine virtual SUSY corrections. For this contribution we find 
relative corrections of the order of 1%, as demonstrated in Fig. 15, which nicely illustrates 
the onset of normal and anomalous thresholds in vertex and box diagrams with rising energy. 
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Figure 15: Corrections from genuine SUSY loops to the cross-section for e+e^ /^r/^r 
relative to the Born cross-section. The corrections are separated into vertex and self-energy 
corrections (a) and box corrections (b). The kinks are generated by normal (single arrows) 
and anomalous thresholds (double arrows). The values are given for the SPSla scenario. 
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Figure 16: Dependence of the relative one-loop corrections (cTq — (7Born)/cBorn to /iR pair 
production on the gaugino parameters Mi, M2 and jj, for ^/s — 500 GeV. The values of the 
other parameters are taken from the SPSla scenario. 
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Figure 17: Electroweak corrections to the cross-section for e~^e~ /^r/^r; relative to 
the improved Born cross-section. Besides the full 0{a) result, contributions from different 
subsets of diagrams are shown. Input parameters taken from the SPSla scenario. 

The variation of the corrections across the [Mi,/i] and [M2,/u] planes has been studied in 
Fig. 16 (a,b). 

All crucial elements have now been collected to present the overall correction of the 
total cross-section for /xr pair production. The parameters of the Snowmass reference point 
SPSla have been adopted again to illustrate the final results. As a function of energy the 
correction normalized to the Born cross-section, defined for the running electromagnetic 
coupling, is displayed in Fig. 17. Moreover, the total correction is broken down to initial- 
state plus final-state QED corrections, the total SM corrections including just lines of SM 
particles in the virtual corrections, and the genuine SUSY corrections introduced above. 
The real photon radiation is treated in a fully inclusive way, i. e. both soft and hard photon 
emission are included. The main contribution to the corrections can trivially be traced back 
to the universal factorizable QED terms which are logarithmically enhanced. However, after 
these dominant effects are subtracted, the remaining QED, weak- loop and genuine SUSY 
contributions still amount to a level of 5%. 

Thus precision measurements of the total cross-sections for slepton pair production, that 
may reach a level of a few per-mille, require the one-loop radiative corrections to be included 
properly. In this way a satisfactory understanding of the expected experimental results will 
be achieved. 
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4.4.2 Selectron production 

In comparison to smuon pair production, the loop calculations to selectron pair production 

are significantly more complex, the reason being twofold. An extra technical challenge is 
introduced by the additional t-channel neutralino exchange mechanisms. These mechanisms 
also give rise to delicate problems for gauge invariant subdivisions of diagram classes and 
the subsequent renormalization procedures. The origin of the problems is the continuous 
flow of charges from the initial to the final states while, at the same time, the Fermi/Bose 
character of the charge line changes in the electron-selectron-neutralino Yukawa vertex — a 
SUSY vertex sui generis. Related problems were encountered first in WW pair production 
via t-channel neutrino exchange as opposed to muon pair production in e'^e~ collisions. 

A transparent example is provided by the process e^e^ ^r^r which is built up solely 
by t-channel neutralino exchanges: 

1. Closed loops of leptons and sleptons implanted in the virtual neutralino lines form a 
gauge invariant subset of diagrams, and so do loops of quarks and squarks. 

2. The diagrams involving massive gauge bosons, Higgs bosons, gauginos and higgsinos 
however cannot be separated from the QED loops in a gauge-invariant manner any- 
more. This follows from a simple argument. The set of photonic corrections to the 
electron-selectron-neutrahno Yukawa vertex is not UV finite, but only so after being 
supplemented by the corresponding virtual photino diagram. Since the photino is not 
a mass eigenstatc, this amplitude is closely linked to the remaining degrees of free- 
dom in the electroweak sector. Thus only the total set of gauge boson/Higgs and 
gaugino/higgsino electroweak diagrams is gauge invariant. 

Nevertheless, as expected on general grounds, soft real-photon radiation regularizes 
the infrared divergences generated by the virtual photon diagrams, and the selectron 
pair cross-section for soft photon radiation factorizes again into the Born term times 
a radiator function. In leading logarithmic order of the soft-photon energy AE, 



d^R^FS?' = d^Born - log (log ^ + 2 log 4 



2 2 

1 + log- 



^ — ^75— log 



(90) 



2/9 1-/9. 

where t and u are the invariant momentum transfers in the t-/u-channel. 

In the same way, in leading logarithmic order, coUinear photon radiation from the initial 
beam fine can be cast in the convoluted form of eqs. (87), (88). Note however that 
virtual initial and final state radiation cannot (even not theoretically) be disentangled 
anymore. 

The infiuence on parameters of the Higgs (and higgsino) sector is rather mild since the 
couplings of these fields to electron-type lines is negligible. Effects of Higgs bosons on the 
self-energies of the Z boson and the neutralinos may naively be expected non-negligible. 
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However, the leading effects of the Higgs boson spectrum on these parameters are only 
proportional to the logarithm of the mass ratios of two Higgs bosons, e.g. oc logMA/Mn. 
As a result, these contributions are naturally suppressed. 

The corrections to the RG improved Born cross-sections for selectron pairs in e'^e" 
and e~e~ collisions are depicted in Fig. 18 (a) and (b). In addition to the full corrections, the 
results are broken down to the individual contributions from closed loops of leptons/sleptons, 
quarks/squarks and the remaining corrections involving gauge bosons, Higgs bosons, gaugi- 
nos and higgsinos, as well as the QED corrections. 

As before, we shall study the influence of the corrections induced through the supersym- 
metry sector at some detail. Even though the genuine SUSY loops are intimately correlated 
with the Standard Model loops, the higher-order effects vary widely over the supersymmetry 
parameter space, measuring the influence of the SUSY degrees of freedom beyond the trivial 
effects due to the masses and couplings of the selectrons produced in the final state. 

A significant influence on the one-loop corrections arises from the electroweak gaugino 
sector, characterized by the parameters Mi, M2 and /i. As an example, the dependence of the 
one-loop corrections relative to the Born cross-section on M2 and /i is shown in Figs. 19 (a) 
and (b). The effects are maximal for small fi due to the higgsino loops affecting the W and 
Z self-energies. [The rapid changes along the diagonal M2 = n are a consequence of the 
level crossings between the Xi states which induces drastic changes in the couplings to the 
electroweak gauge bosons.] 

As outlined earlier, big mass differences between the SUSY sfermions and the corre- 
sponding SM fermions generate large effective Yukawa couplings and thus large superoblique 
corrections to selectron production. This can nicely be illustrated by comparing the squark 
loop effect on the smuon pair cross-sections with the selectron pair cross-section mediated 
by t-channel neutralino diagrams. Beyond the low-mass threshold region, the squark con- 
tributions are rising linearly in the logarithm of the squark masses for selectron production 
while approaching a plateau for smuon production, where gaugino/higgsino lines are absent 
at the Born level, cf. Figs. 20 (a) and (b). 

4.4.3 The identity of SUSY Yukawa and SM gauge couplings 

Selectron-pair production in e'^e~ collisions and, particularly, in e~e~ collisions are excellent 
channels for testing the identity of the SU(2) and U(l) SUSY Yukawa couphngs and the 
corresponding SM gauge couplings. Selectron production in e~e~ coUisions is mediated 
solely by neutralinos, and the same is true for initial-state leptons of equal helicities in e"^e~ 
channels, as evident from Table 1. 

Asymptotically, the t-channel exchange is the dominant mechanism. As a consequence of 
unitarity, the logarithmically leading part of the cross-section for high energies is independent 
of the neutralino mixing parameters after all neutralino exchanges in the t-channel are added 
up: 




(91) 
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Figure 18: Electroweak corrections to the cross-sections (a) for e^e~ ej^ej^ and (b) for 
e~e~ — >• ej^ej^, relative to the improved Born cross-section. Besides the full 0{a) result, 
contributions from different subsets of diagrams are shown. Input parameters taken from 
the SPSla scenario. 
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Figure 19: Dependence of the relative one-loop corrections ((Tq, — (7Born)/cBorn to e^e^ 
production (a) and e^e^ production (b) on the gaugino parameters M2 and /x for y/s — 500 
GeV. The values of the other parameters are taken from the SPSla scenario. 
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Figure 20: Dependence of the one-loop corrections on the soft-breaking squark mass pa- 
rameter Mq [assumed to be universal for all squarks] for (a) fijifiji production and (b) e^e^ 
production. The values of the other parameters are taken from the SPSla scenario, and the 
cms energy is set to y/s = 500 GeV. 
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where g and g' arc the SU(2) and U(l) Yukawa couphngs, respectively. However, very large 
energies indeed would be needed before the asymptotic behavior is reached in practice. 

In this subsection we shall study the sensitivity of selectron production to the Yukawa 
electron-selectron-gaugino couplings and the errors expected experimentally. We assume 
that the masses and mixing parameters of the ncutralinos have been pre-determined in 
chargino/neutralino pair-production, and we properly take into account the expected errors 
in this sector. It is important to note that the mixing parameters affecting the selectron 
production cross-sections depend only on the gaugino/higgsino mass parameters Mi, M2 
and /I. These parameters can be determined in the MSSM from the precision measurements 
of three chargino and neutralino masses, for example, xf, xf and x?) but independently 
of the chargino/neutralino cross-sections [which would re-introduce the Yukawa couplings 
otherwise] . 

To derive the size of the errors on the Yukawa couplings, we perform this study in the 
one-loop approximation. Proper errors will be associated to all terms in the improved Born 
approximation. However, the higher-order corrections can be calculated for the ideal values 
of the parameters in the loops since their errors would affect the errors in the Yukawa 
couplings only to second order which is consistently neglected. Iterative procedures may 
later be employed for the next-level improvements. 

The selectron production cross-sections are computed using the beam polarizations and 
cuts introduced in Section 3.3, and with beamstrahlung and ISR switched on. Besides 
enhancing the statistics, the polarization of both the e~ and e+ beams is essential for iden- 
tifying the chiral quantum numbers of the selectrons [30]. Since only total cross-section 
measurements will be considered here, the one-loop corrections can be included by means of 
a simple K-factor. Besides the neutralino parameters, the cross-sections crucially depend on 
the selectron masses, which can be extracted from threshold scans, cf. Tab. 4. For the neu- 
tralino/chargino masses the following errors are assumed, Sm^o = 50 MeV, m^o = 80 MeV, 
5m-± = 3000 MeV, which are based on a coherent analysis of LHC and LC mass measure- 
ments [55] . From these masses, the gaugino and higgsino mass parameters can be derived, in 
models with two Higgs/higgsino doublets, to which we restrict our general analysis. These 
parameters determine the elements of the chargino and neutralino mixing matrices includ- 
ing estimates of their errors. An error of 1% is assigned to the polarization degree of the 
incoming electron/positron beams. 

As discussed in Section 3, the L- and R-selectron states can be discriminated by con- 
sidering the decay of the selectron into the neutralino X2, followed by the decay chain 
X2 T~^T~ Xi and leading to the final states listed in Tab. 4. It is assumed that each tau 
pair can be identified with an efficiency of = 80%. In addition a global acceptance factor 
of edet = 50% is assigned for potential detector effects, that arc not simulated in this study. 

Taking into account all these statistical errors and systematic uncertainties, the con- 
straints on the SU(2) and U(l) Yukawa couplings of the MSSM, g and g' , are presented in 
Fig. 21 (a) and (b) for the SPSla scenario. The results are based on 500 fb~^ data accumu- 
lated in e+e" collisions at 500 GeV, and 50 fb~^ in e~e~ collisions, respectively. From the 
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(a) e+e-, = 500 GeV, L = 500 fb" 



(b) e-e", y/s = 500 GeV, L = 50 fb"^ 




Figure 21: la bounds on the determination of the supersynimetric U(l ) and SU(2) Yukawa 
couphngs g' and g from selectron cross-section measurements. The two plots compare the 
information obtained from the cross-sections c^ji — (T[e^ e — > Cj^] and cj^l — (^[e^ e — > 
in the e^e~ mode (a) as well as ctrr = a[e~ e~ Cj^] and ai^i^ = a[e~ e~ e'^] 

in the e~e~ mode (b), respectively. Parameters taken from the SPSla scenario. 



overlap regions tlie following expected la errors are obtained: 



e+e" : -^r ~ 0.18%, ^ ^ 1.2%, (93) 

9 9 

e-e- : ^ ^ 0.23%, ? ^ 0.9%. (94) 

9 9 

As evident from these results, the expected sensitivity for the measurement of the Yukawa 
couplings in the e^e~ and e~e~ modes are similar. While e~ polarization is essential for 
disentangling the SU(2) and U(l) couplings, the additional polarization reduces the errors 
by a factor of about 1.4 for a degree of 50%. The result for the determination of the U(l) 
bino Yukawa coupling is comparable to previous studies, Ref. [14,54], while being slightly 
more precise than analyses based on the differential selectron cross-section [16]. 

The precision on the Yukawa couplings expected from selectron pair-production com- 
pares favorably with corresponding measurements in chargino/neutralino pair production, 
for which the errors are of similar though slightly larger magnitude [10,54]. 

Thus it has turned out that one of the basic properties of supersymmetric theories, the 
identity of Yukawa and gauge couplings, can be tested with accuracies down to the per-cent 
and even per-mille level in selectron pair production at e^e~ and e~e~ colliders. 
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5 Conclusions 



We have presented in this report the theoretical basis for high-precision studies of the super- 
symmetric partners to muons and electrons, the scalar smuons and selectrons, at future e+e" 
and e~e~ linear colliders. The theoretical material elaborated in the study is complemented 
by phenomenological analyses of the masses and couplings of these particles. 

Masses: A central target of experiments exploring the properties of super symmetric 
particles is the measurement of their masses. The experimentally observed particle masses 
are connected with mass parameters in the SUSY Lagrangian that encode the breaking of 
supersymmetry and are thus directly related to the basic structure of the supersymmetric 
theory at the TeV scale. Extrapolating these parameters to high scales will allow us to 
reconstruct the fundamental supersymmetric theory. 

Precision measurements can be performed in the clean environment of high-energy lepton 
colliders operating with polarized beams at high luminosity. The masses of scalar sleptons, 
selectrons in particular, can be determined with unparalleled precision from threshold scans 
as the excitation curves rise steeply with the energy above threshold, either with the third 
power in the velocity for smuons or even linearly for selectron channels. 

As a consequence, accurate theoretical predictions for the pair production are required to 
match the expected experimental accuracies. The two theoretical key points in this context, 
non-zero width effects and Coulombic Sommerfeld rescattering effects, have been elaborated 
in detail. Special attention has been paid to preserving gauge invariance in truncated subsets 
of the entire ensemble of Fcynman diagrams describing the final states after the resonance 
decays. The remaining contributions of this ensemble are taken into account as part of the 
SUSY backgrounds, with the additional backgrounds from SM sources added on. 

Based on this procedure, a phenomenological analysis of slepton masses in threshold scans 
was performed, improving significantly on the theoretical reliability compared with earlier 
simulations. While the excitation curves arc characterized by their distinct rise near the 
thresholds, including sub-dominant backgrounds reduces, somewhat, the precision expected 
from previous studies. Nevertheless, a precision of about a few 100 MeV can be expected for 
slepton masses around 200 GeV in general, corresponding to a relative error at the per-cent 
to per-mille level. For the R-selectron mass, an accuracy of even 0.2 x 10~^ can be obtained 
from threshold scans in the e~e~ mode, benefiting from the exceptionally sharp rise of the 
S-wave selectron excitation. 

Moreover, the threshold scans of selectron pair production can also be exploited to extract 
the decay widths to accuracies between 10 and 20%. 

Yukawa couplings: A key character of supersymmetric theories is the identity between 

the Yukawa couplings g{ffV) of the fermions, their supcrpartners and the gauginos, and the 
gauge couphngs g{ffV)/g{ffV) of the fermions and sfermions to the gauge bosons. The 
identity of these couplings is crucial for the natural solution of the fine-tuning problem. It 
must hold not only in theories with exact supersymmetry but also in theories incorporating 
the breaking of supersymmetry, as to ensure the stable extrapolation of the system to energies 
near the Planck scale — one of the defining raisons d'etre for supersymmetry. 
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In passing it may be noticed that the identity of the gauge couphngs themselves in the SM 
and SUSY sectors can be tested at the per-cent level. Smuon pair production is particularly 
suited for extracting the gauge couplings in the SUSY sector as this process is mediated 
solely by s-channel 7 and Z-boson exchanges. 

In contrast, the SUSY Yukawa couplings of the electroweak sector may be probed in 
selectron pair production due to the neutralino t-channel exchange contribution. By carefully 
analyzing statistical errors and systematic uncertainties we could demonstrate that these 
couplings can be extracted from measurements of the total cross-sections in the high-energy 
continuum with a precision of better than the per-cent level. In particular, this slightly 
exceeds the accuracy that can be achieved by other methods, for example in the analysis of 
chargino/neutralino pair production. 

Matching this expected experimental accuracy with its theoretical counterpart requires 
the calculation of the cross-sections to per-cent accuracy. For this purpose the complete 
next-to-leading order one-loop SUSY electroweak radiative corrections were calculated for 
the production of on-shell smuon and selectron pairs. The corrections were found to be 
sizable, being of the order of 5-10%, with genuine SUSY corrections accounting for about 1% 
in pair production, where these corrections can be defined unambiguously and separated 
consistently. 

In all examples analyzed in this report, the production of scalar electrons in e"^e~ an- 
nihilation has been compared to the corresponding processes in e~e~ scattering. The e~e~ 
mode turns out to be particularly favorable for the measurement of the selectron masses 
in threshold scans. In addition it can provide complementary information on the selectron 
Yukawa couplings. 

The detailed analytical expressions for the slepton pair production cross-sections are too 
lengthy to be reported here in writing. The results are implemented in computer programs 
that return the cross-sections for smuons and selectrons to one-loop order for whatever set 
of Lagrangian parameters in the Minimal Supersymmetric Standard Model MSSM is cho- 
sen. The computer codes are available from the web at http://theory.fnal.gov/people/ 
af reitas/. [Technical information on installing and running the programs are given at the 
web site.] 

The theoretical analysis presented in this report for smuons and selectrons is one of 
the cornerstones for precision analyses of the supersymmetry sector in particle physics. 
The precision that is expected to be achieved in future linear collider experiments re- 
quires the analysis of many different channels in parallel — sleptons, charginos/neutralinos 
and squarks/gluinos. The parameters of all these states affect mutually the theoretical pre- 
dictions at the one-loop level so that all the associated production channels must be analyzed 
simultaneously. Such an overall analysis demands complementary and coherent experimental 
action at lepton and proton colliders — a program for which the present analysis is a crucial 
building block^. This ensures finally a self-consistent picture of the SUSY sector at the 
phenomenological level. 

Beyond drain ing a higli-resohition picture of supersymmetry at low energies, the pre- 

^Such a comprehensive study is presently underway, Ref. [56]. 
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cisc determination of these parameters provides the base for exploring the mechanism of 
supersymmetry breaking and the reconstruction of the fundamental supersymmetric theory, 
potentially at scales near the Planck scale. In short, the high precision analyses provide us 
with a telescope for exploring the structure of physics at the scale of ultimate unification of 
genuine particle physics with gravity, as expected to be realized near the Planck scale. 

Appendix 

The theoretical results and the phenomenological analyses presented in this report, have 
been based on the specific reference scenario SPSla for the MSSM, defined in the set of the 
"Snowmass Points and Slopes" [20]. 

The SPSla point is a typical mSUGRA scenario characterized by fairly light sfermion 
masses. If realized in nature, a wealth of experimental information would become available 
on this supersymmetric theory from a linear collider operating in the first phase at energies 
up to about 1 TeV. 

The SUSY parameters of SPSla are defined at the GUT scale for the following universal 
values: 

mo = 100GeV, M1/2 = 250 GeV, ^0 = -100 GeV, tan/3 = 10, > 0. (95) 

The overall set is completed by the Standard Model parameters specified at the electroweak 
scale as 

a(Mz) = 1/127.70, mt = 175 GeV, 

Mz = 91.1875 GeV, = 4.25 GeV, (96) 

Mw = 80.426 GeV, = 1.777 GeV. 

The evohition of the soft SUSY breaking parameters (95) down to the electroweak scale by 
means of the program ISAJET 7.58 leads to the weak-scale parameters listed in Ref. [24]. 
Uncertainties due to the implementation of the renormalization group evolution are not 
relevant for the purpose of the present study, since the MSSM parameters at the weak scale 
are taken as the starting point for our analysis. Using the MSSM soft-breaking parameters 
from Ref. [24] together with the SM parameters (96) the SUSY particles spectrum in Tab. 5 
is obtained. The widths and branching ratios of the particles involved in the decay chains 
of the sleptons can be found in Tab. 2. 

The staus, squarks and Higgs bosons only enter in the loop corrections to smuon and 
selectron pair-production. The masses of the Higgs bosons in the loop contributions must be 
evaluated at tree level, in order to ensure the cancellation of gauge-parameter dependence 
in the next-to- leading order result. On the other hand, some of the background processes 
considered in Section 3.3 involve tree-level Higgs boson exchanges. In order to get rehable 
predictions for these processes, it is necessary to include the large radiative corrections to 
the Higgs masses, which for this purpose were calculated with the program FeynHiggs [57]. 



44 











SUSY particles and 


masses 










142.72 


Mr 


520.50 


ii 


375.74 




96.18 




[89.28] 122.71 


k 


= eL//iL 


202.32 


Ml 


537.20 


h 


585.15 




176.62 




[394.07] 393.56 




= i^e / 


185.99 


(^R 


520.11 


hi 


488.01 


xl 


358.80 


AO 


[393.63] 393.63 




Tl 


132.97 




543.07 


h 


528.23 


xl 


377.87 










206.29 










xf 


176.06 










185.05 










xt 


378.51 







Table 5: Mass spectrum in GeV of SUSY particles relevant for this study for the reference 
scenario SPSla [20]. The Higgs masses are given in Born approximation [square brackets] 
and radiatively corrected. 



Added Note: 

The SPSla set of parameters generates the following predictions for the low-energy precision 
observables: BR[6 57] = 2.7- 10^"^ and A{g^ — 2)/2 = 17- 10^"'^'^. The amount of cold dark 
matter is, with fl^h"^ = 0.18, still compatible with WMAP data but somewhat on the high 
side if they are supplemented by the ACBAR and CBI data. 

Shifting the scalar mass parameter slightly downwards to mo = 70 GeV, but not altering 
any of the other universal parameters in SPSla, drives the value for the density of cold 
dark matter to the central band of the data, Q^h"^ = 0.126, without violating the bounds 
on BR[6 — > sj] and A{g^ — 2)/2 [58]. Slepton, chargino/neutralino masses and branching 
ratios relevant for the present analysis change within so limited a margin that none of the 
conclusions in this report is affected to a significant amount. 
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